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ABSTRACT 

The statistical properties of the ehipticities of galaxy images test the paradigm of 
structure formation and galaxy evolution, and constrain models of galaxy morphology, 
which are key to the removal of the intrinsic alignment contamination of cosmological 
weak lensing surveys. We construct such models based on the halo properties of the 
Millennium Simulation and confront them with a sample of 150, 000 galaxies from the 
COSMOS Survey, covering 4 decades in luminosity and redshifts out to z = 2. The 
ellipticity measurements are corrected for effects of PSF smearing, spurious image dis- 
tortions, and measurement noise. We find that early-type galaxies have a 25 % lower 
intrinsic ellipticity dispersion than late-type galaxies, which is quantitatively repro- 
duced by our best models. None of the samples shows evidence for redshift evolution, 
while the ellipticity dispersion for late-type galaxies scales strongly with absolute mag- 
nitude at the bright end. The low ellipticity dispersions predicted by models based 
on reduced inertia tensors of simulated haloes are generally disfavoured by the ob- 
servations. The fraction of close to circular late-type galaxy images in COSMOS is 
much lower than expected for a sample of circular inclined thick disks, indicating a 
substantial fraction of galaxies with irregular morphology. 
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1 INTRODUCTION 

In the paradigm of hierarchical structure formation tidal 
gravitational torques and shear forces play a central role in 
determining the morphology and angular momenta of dark 
matter haloes over time. These properties affect the way 
galaxies form, evolve and interact with the environment. In 
particular, they strongly impact on the distribution of, as 
well as the correlations between, shapes of the observable, 
luminous parts of galaxies. 

Consequently, the intrinsic shapes and alignments of 
galaxies play a dual role in cosmology: On the one hand 
they constitute a potentially valuable and complementary 
probe of galaxy formation and evolution scenarios, particu- 
larly of the influence of the large-scale gravitational poten- 
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tial in the galaxy's environment. On the other hand, due to 
the scatter and the induced intrinsic correlations, the intrin- 
sic shape properties of galaxies feature prominently in the 
statistical and systematic error budgets of large-scale weak 
gravitational lensing surveys, thereby limiting the accuracy 
obtainable on dark matter, dark energy, or modified gravity 
constraints. 

Both aspects call for a better understanding of the dis- 
tributions and correlations of galaxy shapes, and their de- 
pendence on time, luminosity, environment, merger history, 
and other properties. The large-scale shape corr elations (e.g. 
iMandelbaum et al.l l201ll : IJoachimi et al.l l201ll ). the aUgn- 
ment of satellite galaxies o n small-scales (e.g. iHao et al.l 
l201ll : lHung fc EbelindbOllI'). and the d istribution of galaxy 
ellipticities (e.g. iLeauthaud et al.ll2007r ) have hitherto been 
studied separately. However, a successful model of intrin- 
sic galaxy shape statistics has to explain these observations 
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simultaneously. This paper is the first of a suite in which 
we make a first attempt at constructing such comprehensive 
models and confi-onting them with several sets of new and 
existing observational data, concentrating in this first part 
on one-point statistics of shapes. 

A major goal of this investigation is to establish a 
new, complementary approach to pin down viable models 
of galaxy intrinsic alignments, which, besides reproducing 
the two-point statistics, also have to be capable of predict- 
ing the distribution of galaxy shapes among various galaxy 
populations. Additional constraints would be most valuable 
because current intrinsic alignment constraints are limited 
to z < 0.7, do not extend to deeply non-linear scales (lim- 
ited by galaxy bias measuremen ts), and are still affected 
by large statistical u ncertainties (|Mandelbaum et al.ll201ll : 
Ijoachimi et al.ll201ll . and references therein). 

This is particularly problematic for upcoming weak 
lensing surveys hke KiD^, DE^, LSSlB, and Eucliclil, 
which will have a large fraction if not the bulk of their 
source galaxies at z > 0.7 and retrieve the majority of their 
potentially excellent constraints on cosmology fro m the non- 
linear regime of structu re formation (e.g. Takada fc JainI 
12004 : iLaureiis et al.ll201ll ). Analytical work, N-body simula- 
tions, and observations agree that intrinsic alignments could 
constitute a contaminati on of the order 10 % to weak lensing 
two- point statistics (e.g. Heaven s et al.ll200d : ICatelan et al.l 
l2001j ; Mandelbaum et al. 2006) . 

Hence methods designed to remove or calibrate the in- 
trinsic alignment signal are a necessity but currently have to 
work under minimal assump tions about the form of t he in- 
(se e e.g. iKing fc SchncidCT 2002, 200^; 
\^ I Joachimi fc Schneider 2008. 2009i : 



trinsic correlat i ons 
Bridle fc Kin3 [20m 



BernsteinI l2009l : I Joachimi fc Bridle! I2OI0I : IZhan3 120101 ). A 
robust prediction of the intrinsic alignment contamination, 
which this work aims at, will therefore be a reliable base 
for developing weak lensing survey strategies and tools to 
control intrinsic alignments. 

Much of the early work on the statistical proper- 
ties of galaxy morphologies has focused on inferring the 
three-dimensional shapes of galaxie s from their light dis - 
tributions feinggelilllQSd : iBinnev fc de Vaucouleur3ll98lh . 
iLambas et al.l 1 19921 ) analysed axis ratios of images in the 



APM Bright Galaxy Survey and found significant differences 
in the frequency of small axis ratios between their early- and 
late-type samples. 

The distribution of galaxy ellipticities in fainter samples 
has primarily been investigated to assess statistical error 
limit s on weak lensing m e asure m ents (e.g. Brainerd et al.l 



2007^ . 



1 19961 : iBernstein fc Jarvid I2OO2I : iLeauthaud et al.l 
Nonetheless, some of these results also provided hints at 
clear differences in the ellipticity distributions between dif- 
ferent galaxy populations, and constrained the evolution of 
the dispersion of intrinsic ellipticities with redshift, indicat- 
ing that these measures may add considerable constraining 
power on galaxy shape models. 

Substantial differences in the dispersion of intrinsic el- 
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lipticities between different galaxy populations could have 
interesting implications for measurements of large-scale 
weak gravitational lensing. Forthcoming surveys will cover 
large areas of the sky and will thus be limited by the 
ellipticity noise on medium and small scales from which 
most of the cosmological information is extracted. The 
noise power spectrum is proportional to a^/ng, where a^ 
is the dispersion of the complex ellipticity and n^ the pro- 
j ected number density of galaxi es with shape measurements 
(JBartelmann fc Schneiderll2001^ . 

Therefore certain galaxy samples, appropriately se- 
lected to have low ellipticity dispersion, can beat down sta- 
tistical error limits or become a valid alternative despite 
lower number density. Using such samples might be desir- 
able if e.g. shape measurements became easier, photomet- 
ric redshifts more precise, or intrinsic alignments of galaxy 
shapes either intrinsically weaker or easier to pin down with 
external data. 

In the following we will extract intrinsic ellipticity dis- 
persions and d istributions of ellipt icities from the HST COS- 
MOS Survey (jScoville et al.l l200';i) and confront these mea- 
surements with simulation-based models, with the aim of 
interpreting the statistical properties of galaxy shapes in 
COSMOS, identifying samples that could reduce the noise 
limits of weak lensing surveys, and select realistic models of 
galaxy morphology. In a forthcoming paper we will then use 
the same models to investigate intrinsic shape correlations, 
match them against current observational constraints, and 
predict the intrinsic alignment contamination on planned 
weak lensing surveys. 

As currently it is computationally not yet possible to 
run high-resolution hydrodynamic simulations on a cosmo- 
logical volume, we will rely on a dark matter-only simu- 
lation complemented with a 'semi-analytic' model of the 
galaxy morphology. Our galaxy shape models are based on 
the ha lo properties extract ed from the Millennium Simu- 
lation IjSpringel et al.ll2005l ). which comprises a sufficiently 
large volume to allow for a measurement of large-scale cor- 
relations, but a lso has excellent mass resolution (see e.g. 
iHevmans et al.i i2006 whose simulations have 20 times higher 
particle mass). 

Correlations of dark matter halo ellipticities and an- 
gular momenta among each other and with the large-scale 
matter distribution have been i nvestigated in great de- 
tail with N-body simulations (e.g. Bailin fc Steinmetal2005l : 



lAltav et al.ll2006l : iHahn et al.ll2007i : iLee et al.ll2008l ). We wiU 
supplement this information with multi-band photometry 
and galaxy type classifications from th e semi-analytic mod - 
els of galaxy formation and evolution bv lBower et al.l (|2006l ). 
which enables an accurate selection of galaxy samples for 
comparison with observations. 

As we rely on a dark matter-only simulation, our galaxy 
shape models have to make assumptions about how baryons 
trace the d ark matter. We will follow earlier simula tion- 
based work l|Heavens et al.ll2000l : iHevmans et al. | 20061) and 
analytic intrinsic alignment models (see e.g. ICatelan et al.l 
I2OOII : iHirata fc SeliakI [200J) in assuming that early-type 
galaxies have the same shapes as their dark matter haloes, 
and that late-type galaxies are composed of thick disks per- 
pendicular to the angular momentum of the halo. 

A large number of small-scale, high-resolut i on hydrody- 
namic simulations (|van den Bosch et al.l I2OO2I : ICroft et al.l 
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I2OO9I : iHahn et all I2OI0I : iBett et all I2OI0I : iBettJ |2012| ) have 
been analysed to yield statistical properties of the rela- 
tion between luminous and dark matter, which we incor- 
porate into the models. Moreover the Millennium data in- 
cludes the positions of satellite galaxies, but no shapes 
as the corresponding subhaloes are not sufficiently re- 
solved. Hence we resort to simple models of satellite shapes 
(and alignme nts), partly b ased on the high-r esolution sim- 
ulatio ns by [ Knebe et al.l (i2008); see also J Kuhlen et al.l 
ll200'i1): iPereira et al.l l|2008h : iFaltenbacher et al.1 l|2008l i: 
iKnebe et al.l ( 2010|) for similar investigations into the shapes 
of satellite galaxies and halo substructure. 

This article is organised as follows. In Section[2]we sum- 
marise the main aspects of the underlying simulations and 
the quantities derived therefrom, before detailing in Section 
[3] the modelling of galaxy shapes. We provide an overview 
on the extraction and processing of intrinsic galaxy elliptic- 
ities from the COSMOS Survey in Section H In Section \5\ 
we present the results of our observational analysis and com- 
pare them with various galaxy shape models. We summarise 
and conclude on our findings in Section |6l 

Unless stated otherwise, rest-frame magnitudes are k + 
e-corrected to z = and computed assuming the cosmology 
of the Millennium Simulation (see below) except for a Hub- 
ble constant Hq = 100 /ikm/s/Mpc with h = 1. Magnitudes 
extracted from the Millennium data base are given in the 
Vega system, while all observations use the AB system. If 
direct co mparison is necessary , we resort to the conversion 
tables of iFukugita et al.l (|l996l ). 



2 SIMULATIONS 

2.1 N-Body simulation 

As the basis for our galaxy models we require the shapes 
and angular momenta of the underlying dark matter dis- 
tribu tion, which we obtain from the Millennium Simula- 
tion (|Springel et al.l |2005| ). With a comoving box size of 
500 Mpc/h populated with 2160^ particles of mass rrip — 
8.6 X lO*/i~^M0, the Millennium Simulation provides us 
with a representative sample of the Universe with the reso- 
lution necessary to determine the properties of galaxy-sized 
dark matter haloes accurately. 

The simulation followed the evolution of the matter dis- 
tribution with 64 snapshots from z = 127 to z = using 
the TreePM algorithm of GADGET-2 (|Springelll2005l ) with a 
comoving force softening scale of 5kpc//i. The underlying 
cosmology is a spatially flat ACDM universe with matter 
density parameter f2m = 1 — f^A = 0.25 at redshift zero. The 
z — baryon density parameter is Q^ ~ 0.045, the Hub- 
ble parameter h = 0.73, the power-law index of the initial 
power spectrum rig = 1, and the normalisation of the power 
spectrum as — 0.9. 

These parameters were chosen to be consistent wit h re- 
sults from the 2dF redshift su rvev llPercival et al.ll2002l ) and 
the 1st year data of WMAP fSpergel et al."2003^. More re- 
cent analyses howe ver suggest a significantly smaller value 
of as around 0.8 ()Komatsu et al.l I2OIII : ISchrabback et al.l 
120101 ). The impact of such a change in the normalisation 
of matter fluctuations on the shapes and alignments of 
dark matter halo shapes and angular momenta is not yet 
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Figure 1. Distribution of lialo shapes in the Millennium simula- 
tion. Bottom panel: Shown is the number of haloes as a function 
of their axis ratios (where the axis lengths a ^ b ^ c are the 
square roots of the eigenvalues of the halo inertia tensor) , in the 
redshift ranges < 2 < 0.03 (dotted contours and colour scale) 
and 1.98 < z < 2.14 (solid contours). Contours are plotted at 
1000, 5000, 10000, and 15000 haloes. Note that in this plot pro- 
late haloes reside along the diagonal, oblate haloes along the right 
margin, and spherical haloes in the upper right corner. Top pan- 
els: Probability density of the axis ratios a/c and b/c. 



well understood, but might become particularly relevant on 
small scales where non-linear gravitational physics domi- 
nate. Since in this work we are only interested in the distri- 
bution of galaxy shapes and not their alignments, it should 
be safe to assunre that the high value of erg does not affect 
our results significantly. 

Ray-traci ng through the Mil lennium Simulation was 
performed by iHilbert et al.l ((20091) . We will use those cat- 
alogues which were constructed from 32 light cones with an 
area of 4 x 4 deg^ each (note that we do not require the grav- 
itational shear measurements in this work). This results in 
a mock survey of 512 deg^ out to a redshift of 2: ~ 2.1. Af- 
ter imposing a magnitude limit of F814W < 25 (obtained 
via the semi-analytic models; see Sections 12.31 and 14.21 for 
details), the survey has a mean galaxy number density of 
about 50arcmin~'^. 



2.2 Halo shapes and angular momenta 

We follow iBett et al.l l|2007l ) in identifying bound structures 
in the simulation and in computing their shape and an- 
gular momenta. A dark matter halo is defined as a col- 
lection of self-bound sub-haloes, i.e. single unbound par- 
ticles get discarded. Firstly, groups of simulation parti- 
cles were construc ted with a friends-of-friends algorithm 
(jDavis et al.Ml985l ). followed by the identification of sub- 
haloes a s self-bound structu res within these groups via 
SUBFIND ISpringel et al.ll2005l ). Merger-tree data is then used 
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Figure 2. Minimum number of particles per lialo required for accurate shape measurement. Left panel: Number of halo particles needed 
to achieve a 5 % maximum deviation of the measured axis ratios of the ellipsoid from the input values. The required number is less than 
300 throughout, except for the far top right corner. Right panel: Number of halo particles needed to achieve a 5° maximum deviation of 
the measured orientation of the largest halo eigenvector from the input value. The blue contour indicates the region where the required 
particle number exceeds 300. Note however that for oblate haloes the directions of the two largest eigenvectors become degenerate, so 
that a deviation larger than 5° is tolerable. 



to identify and remove subhaloes that are only transiently 
in proximity to the halo (these are then treated as sepa- 
rate haloes). This procedure removes many of the problems 
associated with friends-of-friends halo identification with- 
out biasing the halo shap e towa rds a spherical boundary, 
as discussed in iBett et al.l (|2007l ) . The resulting halo defi- 
nition corresponds to the 'merger-tree' haloes described by 
iHarker et 



responds t 
aiT (|2006l ) 



The 'shape' of each halo is computed via the quadrupole 
tensor of the mass distribution per unit mass M with com- 
ponents 



M^,^ ^^r 



t,f-l' 1,1^ ! 



(1) 



where A'^p is the number of particles in the halo, and where 
Vi denotes the position vector of particle i with respect to 
the halo centre (defined as the location of the gravitational 
potential minimum) . The eigenvalues and eigenvectors of M 
define an ellipsoid, with the eigenvalues per unit mass giving 
the square semi-axis lengths a^ ^ fe^ ^ c^, and the corre- 
sponding eigenvectors specifying the axis orientations. We 
interpret this ellipsoid as an approximation to the shape of 
the halo. In Fig.[T]we have plotted histograms of the two axis 
ratios of the resulting ellipsoids for redshifts around and 
2. Haloes are preferentially prolate and tend to be closer to 
spheric al at low red shift. 

i Bett et al.l (|2007; ) recommended a minimum particle 
number of A'p = 300 to avoid biases in shape measurement. 
Since this a restrictive condition that would discard more 
than half of the haloes identified in the Millennium Simula- 
tion, we assess the level of accuracy in halo shape measure- 
ment via the following toy model. 

Mock ellipsoidal haloes are generated with a radial 
NFW mass profile and randomly populated with A'p equal- 
mass particles. We assume a concentration of 10 and trun- 
cate the halo at the virial radius. The shapes of these mock 



haloes are then measured via the method outlined above, 
using 100 haloes for each value of A^'p that we test. We vary 
Np from 10 to 10000 in 12 approximately logarithmic steps 
and take input axis ratios in the range 0.1 ^ a/c, b/c ^ 0.9. 

In Fig. [2] we show for every combination of a/c and b/c 
the number of halo particles required to achieve less than 5 % 
deviation of the measured axis ratios from their input values 
and less than 5° deviation of the measured orientation of the 
largest halo eigenvector from the input direction. While gen- 
erally a few tens of particles are sufficient to measure axis 
ratios for triaxial haloes, the requirements on particle num- 
ber become more stringent if two of the axes are of similar 
size, i.e. for strongly prolate and oblate haloes. 

The threshold of 300 particles is only exceeded for close 
to spherical halo shapes, which are not important for our 
analysis as the number of haloes with this shape are small 
(see Fig.[l]), and because these haloes have very low elliptic- 
ity in projection on the sky. The accuracy for the eigenvector 
direction is not met with 300 halo particles or less for oblate 
haloes with a/c > 0.7, but note that in this case the direc- 
tions of the two largest eigenvectors become degenerate, so 
that a deviation larger than 5° is acceptable for our pur- 
poses. The latter result is in good agreement with the limit 
of a/c = . 81 fou nd for haloes with A''p = 1000 particles by 
iBett et all ||2010D . 

The quantity which will eventually be used for further 
analysis is the projected ellipticity of the halo in Cartesian 
coordinates (see Section [3] for details) . Placing the largest 
eigenvector along the line of sight, we compute the accuracy 
in the ellipticity components attainable with 300 halo parti- 
cles. The difference between actual and recovered ellipticity 
varies only weakly with the values of the axis ratios a/c and 
b/c, between around 0.01 for haloes that are strongly ellip- 
tical in projection and around 0.02 for haloes with nearly 
spherical projections. 

When accepting 10 % deviation of the measured axis 
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ratios and 10° deviation of the measured orientation of the 
largest eigenvector, we find that the minimum requirement 
can be relaxed to A^p = 100 with similar accuracy as shown 
in Fig. (2] We will adopt this less stringent limit when mod- 
elling the shapes of central early-type galaxies. 

These lower limits on particle numbers for shape com- 
putati on ar e of the same order as those deduced by Ijind 
l|2002r) and IPereira et al.l l|2008h . the latter paper present- 
ing a similar approach based on mock NFW haloes fo r the 
reduce d inertia tensor M^^"^ — rn-p'Y^^^^ri^fj.ri^v /rf . Ijind 
l|2002f) derive a minimum number of 160 for which ellip- 
ticity correlations are underestimated not more than 5 % 
with respect to high-resolution simulationqj. We agree with 
this work in that sparsely sampled haloes tend to produce 
smaller axis ratios (and hence larger ellipticity on average) 
as well as rapidly increasing uncertainty in the halo orienta- 
tion. The net effect is an underestimation of ellipticity cor- 
relations b y up to a factor of 2 for haloes with 20 particles 
l|jing||2002l ). 

Furthermore we calculate the specific angular momen- 
tum of haloes, 



The merger histories show a clear dichotomoy between 
ellipticals on the one hand and SO and spiral galaxies on 
the other hand, in par ticular with respect to the fraction of 
major merger events (jParrv et al. 2009| ). Since the merger 
history is thought to be decisive for how the morphology of 
the bright part of a galaxy is related to halo properties, we 
hence use the threshold 7?typc = 0.6 to discriminate between 
our early-type and late-type galaxy shape models (see be- 
low). This classification can also be motivated intuitively: 
Lenticular galajdes are disk-dominated systems; their shape 
is thus thought to be determined by angular momentum, 
similar to spiral galaxies. 

The evolutionary models also keep track of whether a 
galaxy is 'central', defined as the galaxy in the most massive 
substructure of a halo at any given time. All other galax- 
ies in the halo are 'satellites', and are treated differently 
with regard to e.g. gas accretion/stripping and orbits (see 
ICole et al.ll200d for details). We adopt this distinction in the 
modelling of galaxy shapes. 



JVp(rvi,) 
i Y^ 



A^p('^vir) ^ 

i— 1 



Vi X Vi 



(2) 



where Vi is the velocity of particle i relative to the halo cen- 
tre of mass velocity. Onl y particles within the virial radius 
are included in the sum. iBett et al.l ([20071) investigate the 
minimum number of particles needed for accurate angular 
momentum calculations by comparing with a low-resolution 
version of the Millennium Simulation. As their Fig. 7 demon- 



strates, below Nr, 



300 the limited resolution of the halo 



causes a sharp upturn in the spin parameter; therefore, we 
adopt this threshold for our computations. 



2.3 Semi-analytic models 

We supplement the information extracted from the simu- 
lation with apparent and rest-frame magnitudes in various 
bands of galaxies hosted by the dark matter haloes using the 
se mi-analytic galaxy e volution model GALFORM in the version 
of Bowe r et al.l (|2006l V Its main updates on previous imple- 
mentations concern the explicit tracking of AGN evolution 
and feedback, the improved modelling of disk instabilities 
and gas cool ing, as well as the us e of the merger tre es by 
iHarker et ahl (2006); for details see lBower et alj (|2006l ') and 

refe rences therei n. 

IParrv et al.l ((20091) classified galaxy morphologies via 
the bulge-to-total ratio of rest-frame 7^-band luminosity, 
Rtype = iff.buige/iK.totai, where the A'-band closely fol- 
lows stellar mass over a wide range of redshifts and is ro- 
bust to uncertainties in modelling details such as redden- 
ing. Defining spiral galaxies via -Rtypo < 0.4, SO galaxies via 
0.4 < -Rtvpe < 0.6, and elliptical galaxies via -Rtypo > 0.6, the 
I Bower et al.l ([2006) models yield a distribution of morpholo- 
gies at low redshift that is consistent with observations. 



3 GALAXY SHAPE MODELLING 

Our modelling of ga laxy shapes adopts the scheme of 
iHevmans et al.l (|200y) in dividing a galaxy sample into late 
types whose shapes are determined by the angular momen- 
tum of the underlying dark matter halo and early types 
whose shapes follow the shape of their haloes. The shapes 
based on halo properties are assigned to the galaxies identi- 
fied by the semi-analytic models as central to the halo, while 
we sample satellite shapes and orientations from distribu- 
tions extracted from the Millennium and other simulation 
works. An overview on the different models presented and 
explored in the following is given in Table [l] 

Note that for this work the orientations of galaxies are 
unimportant because, while orientation has a strong impact 
on the projected shape of an individual galaxy, it is irrele- 
vant for the probability distribution of polarisations of an 
ensemble of galaxies in an isotropic universe. 



3.1 Early- type galaxies 

All central galaxies with -Rtype S5 0.6 and in haloes with 
more than 100 particles are assumed to have the same three- 
dimensional shape as their host haloes. More precisely, we 
project the ellipsoid defined by the eigenvectors and eigen- 
values of the halo inertia tensor onto the plane of the sky 
and treat the resulting ellipse as the shape of the galaxy. 
Let the three unit eigenvectors of the halo inertia ten- 



sor be denoted as Si_i_ — \s 



a:,M) Sy,M) '5||,m} and the absolutc 
values of the semi-axes as o;^ for /i = 1, 2, 3. Then the pro- 
jected ellipse is given by all points x in the plane of the sky 
which fulfil a;^W~^a; = 1, where we have defined 



W 



E 



'X,M^±,/j 



kW 



(3) 



^ Note however that even with several hundreds of particles halo 
shape measurements can still be afflicted with resolution issues, 
e.g. due to unresolved substructure that affe cts the definition of 
subhaloes via SUBFIND jSchneider et aLllJOllh . 



using 






and a = 



3 

E 






(4) 
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Table 1. Overview on models for galaxy shapes. The distinction between central and satellite galaxy is adopted from the semi-analytic 
model. 'Early-type' galaxies have iJtype > 0.6, 'late-type' galaxies -Rtypc < 0.6. The rightmost column contains the identifiers used to 
construct the names of shape models. Note that low-mass galaxies with too few particles in their haloes to make accurate shape and 
angular momentum measurements are assigned random orientations, but otherwise follow the model assumed for the respective central 
galaxy type. 



halo type galaxy type model 



identifier 



early-type same shape as halo; simple inertia tensor 

central " same shape as halo; reduced inertia tensor 

late-type thick disk ± angular momentum; r-odgo-on = 0.25 

" thick disk ± angular momentum; rodgo— on = 0.1 



Est 
Ert 
Sma 
Sth 



satellite 



early-type 



late- type 



shape sampled from MS halo distribution; simple inertia tensor est 

shape sampled from MS halo distribution; reduced inertia tensor ert 

JCnebe_et_al^ ^2008) shape modifications ekn 

thick disk, rpH™_nn = 0.25 sma 



thick disk, r^^„ 



: 0.1 



sth 



iY corresponds to the eigenvector we assume that the statistical halo shape properties of galax- 



Here, sx,m = {sx,^l 

projected along the line of sight. A detailed derivation of 

the foregoing equations is provided in Appendix El 

The galaxy ellipticity i s then defined in terms of the 
complex polarisation e (see iBartelmann fc Schneideii 1200 ll 
for details and other ellipticity definitions), computed from 
the symmetric tensor W via 



ei 

62 



Wu - W22 
Wii + W22 

2Wi2 
Wii + W22 



(5) 



The projection implicitly assumes that the three- 
dimensional light distribution is uniform with a sharp 
cut-off at the edges. We refrain from using more compli- 
cated schemes involving a realistic radial light distribution 
as this could imply variable ellipticity as a function of 
radius, e.g. manifested as isophote twisting. 

We decide to use the shape of the full halo to model 
the galaxy because haloes by definition should be virialised 
and thus would ideally have well-defined, stable shapes. Sub- 
structures, including the most massive ones that host central 
galaxies, are gravitationally bound but not required to be in 
equilibrium, and thus do not necessarily have as well-defined 
boundaries or as stable shapes. 

Many authors studying the morphology of dark mat- 
ter haloes employ the reduced inertia tensor to determine 
shapes, arguing that giving more weight to the inner part 
of a halo reduces the influence of the distribution of sub- 
haloes in the outskirts and produces a better approximation 
of t he shape of the galaxy residing close to the halo cen- 
tre. iBettI (|2012r ) studied the impact of different halo shape 
measurement algorithms using the same data set. Figure 3 
of that work demonstrates that switching from the simple 
inertia tensor (see Equation [l]), to the reduced one increases 
the minor-axis to major-axis ratio by about 25%, with only 
a weak dependence on halo mass. Assuming that a similar 
modification also occurs for the intermediate-axis to major- 
axis ratio, we rescale all semi-axes accordingly to obtain a 
model based on the more spherical haloes resulting from 
reduced inertia tensor measurements. 

For galaxies classified as early type whose haloes have 
Np < 100 we cannot reliably measure their shapes. Instead, 



ies with A^p < 100 are the same as those of more massive 
galaxies. In each redshift slice we construct two-dimensional 
histograms of halo axis ratios from haloes with A'^p ^ 300 
like those shown in Fig.[T] The low-mass galaxies at the same 
redshift are then assigned halo shapes which are randomly 
sampled from these histograms. 



3.2 Late-type galaxies 

All central galaxies with Np ^ 300 and Rtypa < 0.6 are mod- 
elled as circular thick disks whose orientation is determined 
by the angular momentum of the underlying halo. If the ro- 
tation axis of the disk is perfectly aligned with the angular 
momentum vector L — lLx,Ly,L\A , the polarisation of 
the galaxy image is given by 



ei 



62 



cos(26i) 



1 + r^ 
1 „2 



(6) 



sin(26') 



l + r2 ' 
where the polar angle of the image ellipse is computed via 

6 = h arctan ( y^ ) . 

The axis ratio of the ellipse is readily calculated as 



(7) 



ILi, 



~r ^cdgc- 




(8) 



where redge-on is the ratio of disk thickness to disk diameter, 
i.e. approximately the axis ratio for a galaxy viewed edge- 
on. We again assume a uniform light distribution in the disk 
with a sharp cut-off at the perimeter. Moreover we neglect 
any small deviations of the image from an elliptical shape 
in the projection. 

For a disk similar to the one of the Milky Way one 
expects T-cdgc-on to be of the order 0.1, but a representa- 
tive sample of late- type galaxies vie wed edge-on shou l d have 
significant contributions by a bulge. iBailin fc HarrisI (120081 ) 
plot isophotal axis ratios of a large sample of SDSS galaxies 
in their Fig. 4, finding that the smallest ratios for late-type 
galaxies are indeed close to 0.1. Furthermore the distribution 
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quickly drops off below redge-on ~ 0.25, which we therefore 
choose as an alternative value to explore. 

Note that if we had incorporated a spheroidal compo- 
nent explicitly into our models, we would again be faced 
with a radial ellipticity gradient across the projected galaxy 
images, which is beyond the scope of this work. The im- 
pact of bulges implies a distribution of isophotal axis ra- 
tios for galaxies viewed edge-on, where typic al values should 
be brackete d by our two choices of redge-on- iHevmans et al.l 
I1200J, 120061 ) used a similar prescription for late- type galaxy 
models. Disk thickness is accounted for by rescaling ellip- 
ticities as egai = 0.73 ethin disk, which in the edge-on limit 
corresponds to an axis ratio of 0.16, hence lying in-between 
the models we consider. 

Analogous to low-mass early-type galaxies, late-type 
galaxies with A'^p < 300 that have no angular momentum 
information are modelled as randomly oriented thick disks, 
where redge-on has the same value as the model used for 
the corresponding model of central late-type galaxies with 



3.3 Satellite galaxies 

For galaxies residing in the substructures of haloes we do not 
have information about the properties of their dark matter 
distribution. Therefore we have to rely ab initio on assump- 
tions for both the shapes and orientations of satellite galax- 
ies. 

For early-type satellites we proceed in analogy to low- 
mass central galaxies and sample the axis ratios of three- 
dimensional ellipsoids from the histograms obtained for mas- 
sive haloes with shape information at each redshift slice. 
Optionally these axis ratios are rescaled to mimic the use of 
the reduced inertia tensor. The ellipsoids are then oriented 
to point their major axis towards the central galaxy of the 
halo and subsequently projected along the line of sight using 
Eqs. @ to © to yield image polarisations. 

As an alternative model we implement the modifica- 
tions of shapes a nd orientations of subhaloes found by 
iKnebe et al.l (120081 ) in high-resolution dark matter-only sim- 
ulations (KnebeOS model hereafter). Their measures of 
triaxiality and sphericity of the satellite population can 
be converted to average axis ratios {a/c}^^^ = 0.80 and 
{blc)„,t = . 90, w hich is in agreement with earlier works. 
iKnebe et al.l (|2008l ) calculate the corresponding mean quan- 
tities for central galaxies from previous publications, obtain- 
ing (a/c)^^^ = 0.66 and (b/c)^^^ = 0.76, i.e. the satellite 
galaxies have more spherica l shapes than central galaxies 
(see also lKuhlen et al.ll2007l ). 

We account for this by rescaling all axis ratios {a-fcj^cn 
that are sampled from the histograms of central halo shapes 
via 



(a/c) = 1 - ^ {"{""j^-^ [1 - (a/c) 1 



(9) 



and likewise for b/c. Note that this formula is applied to 
res ults for the r educed inertia tensor as this was also used 
bv lKnebeetHI |2008). 

Late-type satellites are assumed to be thick circular 
disks (with their angular momentum perpendicular to the 
line connecting the position of the satellite with the cen- 
tre of the halo) with the same properties as central disk 
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Figure 3. Histograms of COSMOS LRG (red lines; multiplied 
by a factor of 10 for easier inspection) and non-LRG (black lines) 
galaxies as a function of F81AW apparent magnitude (bottom 
left), redshift z (top left), and rest-frame magnitude My (top 
right). 



galaxies, i.e. we create two models with rcdgc-on = 0.1 and 
rcdgo-on = 0.25, respectively. 



4 COSMOS GALAXY SHAPES 

4.1 Data 

We base our anal ysis on the HST COSMOS Survey 
(jScoville et al.l 120071 ) which is the largest space-based sur- 
vey to date and comes with excel lent photometric red- 
shift information (|llbert et al.l 120091 ) . To estimate the in- 
trinsic polarisation dispersion of galaxies as a function of 
brightness, galaxy type, and redshift, we make use of the 
shape catalogue produced for the weak lensing analysis by 
Schrabback et al. (2010). 

We use only those galaxies for which the absolute V- 
band magnitude (estimated only for the COSMOS-30 sub- 
sample) and the photometric redshift are available. For those 
galaxies the characteristic size in terms of the half-light ra- 
dius, the observed magnitude in the _F814W^ band including 
an error estimate, and a LRG flag are also given. The lat- 
ter identifies galaxies with Mv < — 19 and a photometric 
type classifying them a s 'elliptical' (including SO) as LRGs 
IjSchrabback et al.ll2010l ) . To allow for an accurate treatment 
of noise, we will also impose a magnitude cut _F814W^ < 25, 
which provides us with a total number of 150, 000 galaxies 
of which about 5, 500 are LRGs; see the histograms in Fig.|31 

Our COSMOS catalogue is incomplete for very ex- 
tended, and hence bright, objects because galaxies which do 
not fit well into the postage stamps used for shape measure- 
ment are discarded. Galaxy images with a half-light radius 
of Th = 0.75" or larger were excluded before analysis, and a 
substantial fraction of objects with slightly smaller r^ were 
subsequently flagged as having problematic shape measure- 
ments. This removes in particular bright galaxies, including 
LRGs, at low redshifllj. However, unless this population has 



As these objects have large angular size, PSF effects are likely 
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considerably different ellipticities than less extended galax- 
ies with similar redshifts and luminosities, we do not expect 
a significant impact on our results. 

Gala xy ellipticities we re determined with the KSB+ 
pipeline (jErben et al.ll2001^ ■ which measures weighted sec- 
ond brightness moments of galaxy images. Great care has 
been taken to remove the effects of PSF smearing as well 
as spurious ellipticities due to image distortions, including 
spatial and temporal variations thereof. A detailed account 
of the shape m easurement and te s ts of various systematics 
can be found in lSchrabback et al.l l|2010f ). 

The KSB scheme is designed to provide accurate and 
high signal-to-noise estimates of the gravitational shear 
rather than measure intrinsic galaxy ellipticities, giving a 
strong weight to the inner parts of the galaxy image. Note 
again that real galaxy images generally do not have a sin- 
gle ellipticity, but both absolute value and position angle 
of the ellipticity can be a function of radius. However, we 
will work with the single weighted polarisation estimate, 
once appropriately corrected, that is provided in the COS- 
MOS shear catalogues. This needs to be kept in mind when 
confronting the COSMOS measurements with our simple, 
single-ellipticity galaxy shape models. 

KSB produces an estimate of the galaxy ellipticity e, 
whic h is readily converted to polari sation vi4j e = 2e/(l -|- 
\t\^) i|Bartelmann fc Schneid er'2001'). While both definitions 
of galaxy ellipticity are equivalent, we choose c in this pa- 
per as the measurement noise correction (see below) be- 
comes slightly more convenient computationally, and as the 
ellipticity distributions display their features more clearly in 
terms of e. 

The photometric redshift quality is excellent with a^ ~ 
0.01(1 + z) for galaxies with Subaru i^ < 24 and z < 1.25, 
degrading to a^ ~ 0.06(1 + z) for the fainter galaxies 
at the majcimum redshifts in our analysis around z — 2 
Ijllbert et al.ll2009l V Photometric redshift scatter leads to the 
smoothing of any features in the signal when considered as 
a function of redshift or rest-frame magnitude (if the pho- 
tometric redshift estimate is used to compute the distance 
modulus). However, as will be shown below, the modelled 
signals do not show strong features and are smoother than 
the observed ones even without including the effect of pho- 
tometric redshift scatter, so that its influence can safely be 
neglected. 



4.2 Method 

The polarisation dispersion is computed as 



o-e = y/{ee*) = 



\ 



N 



(10) 



where A'^ is the number of galaxies in a given bin. We choose 
this quantity as our default measure because it provides 
information about the distribution of e in a compact way 



to be negligible, so that the shapes of these galaxies could readily 
be obtained from shallower surveys with less favourable image 
quality. 

Note that this relation holds for the ellipticities defined with 
unweighted (or appropriately corrected) brightness moments. 
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Figure 4. Polarisation dispersion CTj as a function of apparent 
F81AW magnitude. Top panel: The blue dot-dot-dashed line cor- 
responds to (Te measured from the 'raw' polarisations in the shear 
catalogue, the red dashed line to Ue obtained from the KSB shear 
estimates, and the black solid line to the noise-corrected tJe. We 
also process the raw polarisations directly, resulting in a noise- 
corrected ae shown as violet dot-dashed line, and Ce that has sub- 
sequently been corrected for the circularisation due to the Gaus- 
sian weight in the brightness moments, see the black dotted line. 
Note the good agreement between the two black curves at bright 
magnitudes. Bottom panel: Dispersion CTe of th e measurement 
noise, as estimated from the images according to lHoekstra et al.l 
£200Cj), see the black dotted line, and from our Fisher matrix for- 
malism, see the grey solid line. The two approaches show very 
good agreement. 



and, besides, is of relevance to the noise computation in 
weak lensing two-poi nt statistics (e.g. lSchneider et al.ll2002l : 
I Joachimi etallbOOSl ). 

The polarisation dispersion measured from COSMOS 
data is composed of the intrinsic polarisation dispersion and 
contributions from measurement noise. To correct for the 
latter, we use a m odified version of the F isher matrix ap- 
proach proposed bv lLeauthaud et al.l ([20071). Assuming that 
the light distribution in the image can be described by a 
bivariate Gaussian whose covariance is given by the second- 
order brightness tensor, one obtains an estimate of the error 
on the brightness moments which can be propagated into 
an error on ae. Since the Fisher matrix provides us with an 
expectation value, we do not need to revert to the actual 
galaxy images but require only the brightness moments, the 
apparent magnitude, and the S/N for each galaxy as input. 
Full details of this procedure are provided in Appendix IB II 

In a subset of the COSMOS field we also calculate the 
measurement noise directly from the background noise of in- 
dividual galaxy images, using th e method outlined in the ap- 
pendix of lHoekstra et al.l(|2000r ). The resulting measurement 
noise contribution to Ue at different F814:W magnitudes is 
shown in the bottom panel of Fig.U] together with the pre- 
diction of the Fisher matrix formalism. The agreement of the 
two approaches is very good down to our magnitude limit, 
which is remarkable given that our version of the Fisher 
matrix calculat ions does not have any free parameters, as 
opposed to the iLeauthaud et al.l ((20071) formalism. 
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In Fig.m top panel, we have plotted Ue as determined 
from the output of the KSB shear estimation pipeline, as 
well as after subtracting in quadrature the measurement 
noise from the Fisher estimates. As expected, the difference 
between the two curves is close to zero at bright magnitudes 
where galaxies have high signal-to-noise wherefore measure- 
ment noise is negligible. This difference gradually increases 
for fainter magnitudes until the dispersion of the measure- 
ment noise reaches about 0.2, which marks a 10 % contribu- 
tion to o-^ at F814W = 25. 

W e compare this result to Fig. 17 of iLeauthaud et al.l 
l|2007h who used galaxy shape measurements of COSMOS 
jalajci es based on the method proposed by [Rhodes et al.l 
20001 ). Note however that the authors plotted the mean 
dispersion of the ellipticity components ei,2 rather than the 
dispersion of the total polarisation e as done in this work. 
Generally, we find good agreement but obtain overall slightly 
lower values of the noise-corrected ellipticity dispersion, de- 
spite smaller noise-corrections (for _F814W^ < 21.5 we find a 
sh ear dispersion per com ponent of 0.23 compared to ~ 0.24 
in lLeauthaud et al.ir2007l ). 

Small differences in the ellipticity dispersion are not 
unexpected due to the different shape measurement meth- 
ods used. For instance , in the KSB implementation of 
ISchrabback et al.l (J2010l ) the shear tensor is individually de- 
termine d for each galaxy instead of for ensembles of galax- 



ies as m 



Leauthaud et al.l ([20071). Very similar to our Fig.|4l 

ILeauthaud et al.l (20071) also find a slight increase of ae as 



a function of apparent magnitude, even after noise correc- 
tion, which, if physical, could be caused by changes in the 
galaxy population (e.g. due to redshift evolution or changes 
in the fraction of early- and late-type galaxies) at these mag- 
nitudes. 

To cross-check the polarisation measurement by means 
of KSB shear estimates, we also determine a^ via an alter- 
native route: We start from the 'raw' polarisations directly 
obtained from the observed brightness moments and cor- 
rect them for measurement noise via the methods outlined 
above. The resulting polarisations are still affected by PSF 
smearing as well as spurious ellipticity introduced by image 
distortions. One of the main advantages of employing shear 
estimators to measure intrinsic galaxy shapes is that these 
PSF effects have already been removed, but we can proceed 
with the alternative method for bright and extended objects 
whose apparent size is large compared to the PSF width. 

Moreover, the brightness moments are computed with a 
circular Gaussian kernel (see ISchrabback et al.l l20ld ) , with 
the consequence that, even for bright galaxies, the measured 
absolute value of the polarisation |e| will be significantly 
smaller than the 'true' value, as e.g. measured from isophotes 
(which is impossible with sufficient accuracy for the majority 
of faint and small galaxies in COSMOS). In Appendix lB2l we 
show how to correct for the circularisation of the brightness 
moments by quantifying the effect using analytic Sersic light 
profile^ 



* Note that in the KSB scheme this step corresponds to the di- 
vision by the 'shear tensor' which accounts for the circularisa- 
tion due to the Gaussian kern el as well as due to the PSF (e.g. 
iBartelmann fc SchneideJbOOll) . 



The dotted lines in the top panel of Fig. |4] illustrate the 
steps outlined above. The raw and noise-corrected polari- 
sations yield very small CTc < 0.4, which is caused by the 
circularisation effects of the PSF and the circular Gaussian 
weight in the brightness moments. Correcting for the latter 
indeed increases ae to be in very good agreement with the re- 
sult of our default approach. Beyond F814W ~ 22 however, 
the curves diverge quickly as the PSF smearing for these 
fainter and smaller galaxy images becomes an increasingly 
important effect. 

It is also instructive to compare our results for the raw 
polarisation with those of Fig. 20 in iHoekstra et al.l (120001 ) 
who found that the dispersion is constant with apparent 
magnitude, which they traced back to the fact that increas- 
ing PSF effects at fainter magnitudes happen to be exactly 
balanced by an increase in measurement noise. Their mea- 
surements were based on shallower HST imaging, so that for 
COSMOS we expect the PSF effects to wi n over noise for 
F814W > 22 where the noise estimates by iHoekstra et al.l 
(2000) become substantial. This is indeed the case as evident 
by the downturn of ae at FSIAW ^ 22. 

Once we have determined the final, noise-corrected po- 
larisations for all galaxies, the polarisation dispersion is cal- 
culated as a function of rest-frame magnitude Mv and red- 
shift z. Additionally, we divide the galaxies into an LRG and 
a non-LRG sample because, given the different morphology, 
we expect a different a^ for the two samples. 

The error on the dispersion is estimated by bootstrap- 
ping from 50 catalogues containing the same number of 
galaxies as the input catalogue (after the selection criteria in 
magnitude and morphology have been applied). The corre- 
sponding statistics from the simulations are constructed as 
follows. Galaxies are randomly resampled from the simula- 
tion catalogues according to histograms in the {Mv,z) grid 
for the COSMOS LRG and non-LRG samples, respectively. 
While Mv rest-frame magnitudes are directly available from 
the semi-analytic models, we choose to reproduce the LRG 
criterion by selecting all galaxies with -Rtype > 0.4, i.e. lentic- 
ular and elliptical galaxies according to the classification by 
iParrv et aj] (|2009| ). in addition to the cut My < -19. 

We build 10 catalogues for each of the 32 lines of sight 
with numbers of galaxies chosen such that the total in each 
line of sight is of the same order as the number of galaxies in 
the COSMOS catalogues. This ensures that the statistical 
constraints from the simulations are at least as good as from 
the observations, and that repeated draws of galaxies from 
regions in the {Mv, z) grid with sparse sampling are kept to 
a minimum. We show and analyse only bins which contain 
30 galaxies or more after applying all cuts. 

In addition to the selection in Mv and z we have to ac- 
count for the apparent magnitude cut in ^814^. As this fil- 
ter is not available in our simulation catalogues, we approx- 
imate it by the SDSS bands i and z which ro ughly cover the 
same wavelengths. Resorting to the tables of lFukugita et al.l 
(|l996l ), we indeed find that the colour (i + z)/2 - F81AW 
evolves little with redshift. It is essentially constant at 0.31; 
only for early- type galaxies above 2 = 0.5 does this value 
start to decrease moderately. With this conversion we im- 
pose the cut F81AW < 25 on our models. We also test a 
brighter magnitude limit F814W^ = 24, but find that a^ 
remains basically unchanged for all models. 
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Figure 5. Polarisation dispersion CTe as a function of redshift (left panels) and rest-frame V-band magnitude (right panels). The top 
panels display results for the non-LRG sample; the bottom panels for the LRG sample. The dispersion based on the noise-corrected 
polarisations from COSMOS are shown as black filled squares for an apparent magnitude cut at i<'814H^ = 25, and as black open squares 
for a cut at FUliW = 24 (not shown in the bottom panels as almost coinciding with the solid squares). The corresponding results for 
the simulation-based models are shown as red circles (brown diamonds) when basing early-type models on the simple inertia tensor 
[Est;est] and using late-type galaxies with r^dge— on = 0-1 [Sthjsth] (rcdgc— on = 0.25 [Sma ; sma] ) . The same late-typo models combined 
with elliptical galaxies based on reduced inertia tensor measurements including the KnebeOS model for satellites [Ert ; ekn] are given by 
the blue downward triangles (?'g;jgg_ou = 0.1) and the violet diamonds (r^dge-on = 0.25), respectively. Discarding the modifications due 
to the KnebeOS model changes the model from the violet diamonds to the green upward triangles. Note that the LRG sample has been 
defined with a magnitude cut My < —19, as indicated by the black dotted line. Error bars are shown throughout, but remain much 
smaller than the size of the symbols. 



5 RESULTS 

5.1 Polarisation dispersion 

The resulting polarisation dispersion ae for the COSMOS 
data as well as for various galaxy shape models based on 
the Millennium Simulation are shown in Fig.[Sl as a func- 
tion of l/-band rest-frame magnitude and redshift. We dis- 
tinguish between models in which the shapes of early-type 
central and satellite galaxies are either computed from the 
simple [Est ; est] or the reduced [Ert ; ert] inertia tensor (see 
Table [1] for a list of the model identifiers). Moreover the 
KnebeOS model for early-type satellites is used in combi- 
nation with reduced inertia tensor shapes for central early 
types [Ert ; ekn] . For central and satellite late-type galax- 
ies we vary the edge-on axis ratio between rcdgo-on = 0.1 
[sthjsth] and rcdgo-on = 0.25 [Sma; sma]. 

The most striking feature of the observational results 
is a clear difference between a^ for the LRG and non-LRG 
samples, the latter being significantly higher for all magni- 
tudes and redshifts considered. All simulation-based models 
reproduce this dichotomy qualitatively, but only the model 
Est-Sth-est-sth (red circles) agrees very well with the 
COSMOS data in overall amplitude for both LRG and non- 
LRG samples for the redshift as well as the luminosity de- 
pendence. 

Models of the LRG sample are most sensitive to the as- 



sumptions of early-type galaxy shape modelling, while the 
(7e amplitude for non-LRG models is mostly determined by 
late-type modelling details. This is not entirely obvious as 
the LRG models include SO galaxies which follow our late- 
type disk models, whereas the non-LRG models include sim- 
ulated early-type galaxies at My > —19. Incorporating the 
additional rounding of early-type satellites as suggested by 
the KnebeOS model only marginally lowers CTe at low redshift 
and faint My where a larger fraction of satellite galaxies 
contributes to the signal. 

All models feature relatively little variation of (Te with 
redshift and luminosity. Consequently, even the best-faring 
model deviates weakly but still significantly from the COS- 
MOS data in several bins. The most apparent discrepancy 
is observed for the My dependence of the non-LRG sample, 
where the observations favour a low value of a^ for the most 
luminous galaxies, a pronounced increase in the dispersion 
by about 0.15 towards fainter magnitudes, and a subsequent 
decrease in a^ for My > —19, whereas the models remain 
close to constant. 

It is interesting to note that the considerable shift of 
axis ratios of massive haloes in the Millennium Simulation 
towards rounder shapes at low redshift, as shown in Fig.[T] 
is not reflected in the redshift dependence of a^. On the 
contrary, the majority of models suggest a weak trend to- 
wards lower (Te at high redshift for the LRG sample which 
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Figure 6. Polarisation dispersion (Je as a function of rest-frame 
V-band magnitude and redshift for the COSMOS non-LRG sam- 
ple with FS14W < 25. There is a clear tendency towards lower o-f, 
for luminous galaxies with My J; —22 in a broad redshift range 
0-5 < 2 < 1-8. Note that a cut at FSMW = 24 preferentially 
depopulates the bins at lower luminosity and hence enhances the 
impact of the low-cre region at a given redshift. 



is dominated by early- type galaxy models. This could be 
due to our selection criteria, e.g. the bulk of haloes which 
evolve towards spherical shapes a.t z ~ could host late- 
type galaxies whose shapes are not based on the shapes of 
their underlying haloes. 

The non-LRG signal from COSMOS also displays small 
variations in the redshift dependence, with a shallow max- 
imum at z ~ 0.5 and a minimum around z ~ 1.5. Similar 
small-amplitude variations of the redshi ft dependence were 
also observed bv lLeauthaud et al.l ((20071). Changing the ap- 
parent magnitude cut to FSIAW — 24 considerably deep- 
ens the minimum at high redshift and also increases the 
variations in the luminosity dependencqj. Together with the 
two-dimensional distribution of CTe shown in Fig.|6l this al- 
lows for an interpretation of the amplitude variations for the 
COSMOS non-LRG sample. 

While mostly fluctuating around (Je ~ 0.6, the polar- 
isation dispersion strongly decreases down to a^ ~ 0.4 for 
Mv < —21 over almost the complete redshift range that the 
observations cover, i.e. for 0.5 < 2 < 1.8. The brighter mag- 
nitude limit depopulates the bins at the low-luminosity end 
of the plot, thereby rendering the contributions by the region 
of low (Te more important at a given redshift. This feature 
explains the drop at bright My as well as the smaller ae 
in the range 0.5 < z < 2, both becoming more pronounced 
with FSUW < 24, as seen in Fig.El 

It is not unreasonable that very bright disk galaxies, 
which preferentially reside in high-density regions, have dis- 
tinctively different morphological properties compared to 
field galaxies. It is interesting to note that CTe for Mv < —22 
reaches values very close to those predicted by the models 
that include disk galaxies with redge-on = 0.25, which was 



^ Note that this change in limiting magnitude affects the mod- 
els, as well as the COSMOS results for the LRG sample, only 
marginally, so that we do not show the corresponding data in 
Fig.[5] 



motivated by the existence of bulges. Hence the observed de- 
crease in (Te with increasing luminosity could actually be a 
transition from galaxies with no or subdominant bulges (and 
thus rcdgo-on ~ 0.1) to galaxies with a prominent spheroidal 
component, with corresponding links to the merger history 
of these galaxy populations. 

However, we caution that the split into LRG and non- 
LRG sample was based on a cut in the photometr ic type 
between SO and earlier , and Sa and later (Schrabback et al.l 



l20ld : lllbert et al.1l2009l : iPolletta et "ai]|2007^ ■ Hence it is pos- 
sible that the two samples simply become very similar as re- 
gards their Hubble type at the highest luminosities. A finer 
binning in morphological type could test this scenario, which 
we leave to future work. 



5.2 Polarisation distributions 

The limitations of our simple shape models become apparent 
when considering the full distribution of the absolute value 
of 4ij- In Fig. Owe have plotted the probability distributions 
of |e| for the full LRG sample and non-LRG galaxies in the 
magnitude range — 19 < Mv < — 17 (which in this range is 
the full galaxy sample). 

Since the correction for measurement noise outlined in 
Section 14.21 is statistical in nature, we cannot de-noise the 
polarisation measurements of individual galaxies. Thus we 
have plotted the distribution of |e| as produced by the shear 
pipeline, i.e. including all corrections for PSF effects, but 
also measurement errors due to noise. To allow for a fair 
comparison, we devise a simple noise model and apply it to 
the model polarisation distributions as follows. 

We divide the COSMOS galaxy samples used to pro- 
duce Fig. [7] into ten equidistant bins in |e| and calculate for 
each bin the mean dispersion of the measurement error ac- 
cording to our Fisher matrix formalisig^^l. Each polarisation 
from the simulation-based catalogues is then modified by 
a shift randomly drawn from a zero-mean Gaussian with a 
width equal to the corresponding measured mean measure- 
ment error dispersion. Resulting polarisations with \e\ < 
or |e| > 1 are discarded and resampled, so that, effectively, 
the scatter due to measurement noise preferentially increases 
the ellipticity of nearly round objects and decreases the ellip- 
ticity of galaxies viewed nearly edge-on, as one would expect 
in reality. 

The impact of this noise model on the probability distri- 
butions of |e| is illustrated by comparing the two red curves 
of the Est-Sth-est-sth model in Fig.[71 where the dotted 
one is without noise. The scatter due to noise smooths sharp 
features in the distributions, such as the peak due to disk 
galaxies viewed edge-on at |e| < 1 or the strong concentra- 
tion of circular galaxy images in the non-LRG sample at 
|e| =0. 

In the case of the LRG sample the models that were 



^^ Note that CTj, is not the scatter of this distribution but the 
standard deviation of the distribution of the complex polarisation 
e which has zero mean in both its real and imaginary part. The 
distribution shown is still a useful indicator as a^ is computed 
from |e|, see Equation mOI I. 

^^ The mean dispersion of the measurement error generally in- 
creases with |e|, with overall higher values as well as a stronger 
increase for the non-LRG sample. 
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Figure 7. Distribution of the absolute value of the polarisation 
|e|. Top panel: Probability density of |e| for the non-LRG sample 
in the rest-frame magnitude range —19 < My < —17. The curve 
for the COSMOS polarisation as output by the KSB pipeline (i.e. 
without measurement noise corrections) is shown in black. The 
solid and dashed curves resulting from the simulation-based mod- 
els with measurement noise added are assigned the same colour 
coding as used in Fig.fS] The dotted red curve represents the re- 
sult for the Est-Sth-est-sth model without noise. Bottom panel: 
Same as above, but for the full LRG sample. 



the non-LRG sample (see Fig. [5]) may be coincidental, given 
thes e substantial deviat i ons of the |e| distributions. 

IBernstein fc Jarvid ([2003) presented similar intrinsic 
polarisation distributions for bright, low-noise galaxies in 
the CTIO lensing survey. Interestingly, their low surface 
brightness sample, dominated by spiral galaxies, also con- 
tains few near-circular galaxies, with a strong decline of the 
probability distribution of |e| below |e| = 0.1. The authors 
concluded that the disks of these late-type galaxies are not 
perfectly circular, which indeed seems to be the most plau- 
sible explanation. 

Changing the axis ratio of the light distribution of a 
late-type galaxy viewed face-on from 1 to 0.9 already trans- 
lates into a minimum ellipticity of |e| ^ 0.1. This could 
readily be achieved in practice by the presence of luminous 
substructure such as giant star-forming regions or blended 
satellite galaxies. Furthermore it should be kept in mind that 
KSB-like methods put a strong weight on the inner parts of 
a galaxy in measuring brightness moments, so that features 
like prominent bars may play a non-negligible role. Deep 
high-resolution data indeed reveal an increase in the frac- 
tion of asymmetric and irregular galaxies at higher redshift 
^Abraham et al.lll996l : lBundv et al.ll2005l '). so that a lack of 
low- ellipticity galaxy images in COSMOS is not unexpected. 

Both CTIO and COSMOS samples contain a substan- 
tial yet subdominant fraction of early-type galaxies, which 
- if not correctly predicted by the Millennium semi-analytic 
models - could lead to a deviation of the shape models. 
But since both obse rvational data sets (see th e high surface 
brightness sample of lBernstein fc Jarvisll2002l ) as well as our 
best models agree in that early-type galaxies have predom- 
inantly low-ellipticity images, neither a higher nor a lower 
fraction of early-type galaxies in the non-LRG sample could 
render a substantial amount of modelled galaxy shapes with 
lei < 0.3 into ones with lei > 0.5. 



close to the COSMOS results in terms of Ue in Fig. [5] per- 
form reasonably well in reproducing the full distribution of 
|e|. However, the model distributions are more peaked and 
therefore produce slightly smaller values of <je than observed 
with COSMOS. The models that include rounder early-type 
galaxy shapes based on reduced inertia tensor measurement 
for the dark matter haloes are inconsistent with the data 
at high significance, irrespective of what is assumed for disk 
galaxies. 

The |ej distribution of the non-LRG sample peaks at 
|e| ~ 0.8, remaining at high probability up to |e| < 0.95, and 
is hence largely different from the LRG distribution which 
obtains its highest values in the range 0.1 < |e| < 0.4. All 
models are clearly discrepant with the observational data in 
the non-LRG sample, the strongest deviation occurring at 
|e| < 0.2 where the models are dominated by face-on disks, 
whereas the fraction of near-circular galaxy images in this 
COSMOS sample is very small. Even a more sophisticated 
noise model that e.g. does not rely on the assumption of 
Gaussian scatter is unlikely to remedy this difference. The 
Est-Sth-est-sth model fares best by predicting the small- 
est peak at low |e| and a relatively high fraction of highly 
elliptical galaxy shapes, but the good agreement between 
the amplitude of a^ between this model and COSMOS in 



6 CONCLUSIONS 

In this work we analysed the statistical properties of galaxy 
ellipticities, confronting a sample from the HST COSMOS 
Survey with a suite of models based on the dark matter halo 
and galaxy properties provided by the Millennium Simula- 
tion. The galaxy shape models differentiate between central 
and satellite as well as early- and late-type galaxies, and 
incorporate additional information on the link between lu- 
minous and dark matter from other simulations. We con- 
firm earlier work in that at least 300 particles per halo are 
required to measure accurat e three-dimensio nal shapes (as 
well as angular momenta, see lBett et al.ll2007l ). implying de- 
viations in the Cartesian components of the projected ellip- 
ticity of at most 0.02. 

Intrinsic galaxy ellipticities for a sample of 150, 000 
COSMOS galaxies were extracted from gravitational shear 
estimates based on second brightness moments, taking full 
advanta ge of the scrutiny and systematics testing under- 
taken bv lSchrabback et al.l (|2010l ). particularly with respect 
to the smearing of the PSF and spurious ellipticities intro- 
duced by telescope and camera. We demonstrated that the 
resulting polarisation dispersion ae is robust with respect to 
the methods used to correct for the effects of the circular 
kernel in the brightness moments, as well as measurement 
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noise. To estimate the latter, we devised a Fisher matrix 
formahsm that yields accurate results without the need to 
resort to observed or simulated images. 

Splitting the COSMOS galaxies into an LRG and a non- 
LRG sample, we detect a significant dichotomy in the po- 
larisation dispersion between early- and late-type galaxies, 
the latter having a^ that is 0.1 to 0.15 larger on average. A 



notable exception are the brightest (rest-frame My 



< 



-21) 



non-LRG galaxies which feature low Ue close to the values 
found for early types, over a broad range in redshifts out 
to z ~ 1.8. This effect could hint at an environment depen- 
dence, with lower ellipticities found in high-density regions. 
In the case of disk-dominated galaxies, the decrease in Ue 
can be explained by the presence of more prominent bulges, 
reducing the number of (edge-on) galaxies with high ellip- 
ticity close to unity. 

A caveat is that the decrease in ae for bright late-type 
galaxies might also be driven by a transition in galaxy type 
within the broadly defined non-LRG sample towards ear- 
lier types at high luminosity. A future investigation which 
splits the galaxies into finer galaxy type bins and addition- 
ally measures a dependence of CTe on local density will be able 
to distinguish between these two proposed explanations. 

Apart from effects caused by the population of bright 
low-ellipticity late-type galaxies, we find no evidence for a 
redshift evolution of a^ in any sam ple, which is in agreement 
with earlier work (jLeauthaud et al. 2007 ). Studying the dis- 
tributions of absolute values of the polarisation, we identify a 
strikingly low fraction of close to circular galaxies (|e| < 0.2) 
in the non-LRG sample. We agree with iBernstein fc JarvisI 
l|2002r ). who made a similar observation in a brighter galaxy 
sample, that the most likely explanation is a non-circular 
light distribution of disk galaxies viewed face-on, caused e.g. 
by luminous substructures, plus a contribution by irregular 
galaxies. 

The considerable variation of CTe with galaxy properties 
(up to ~ 0.3 in extreme cases; see Fig.O suggests that a 
better understanding of the physics that drive the observable 
shapes of galaxies might help optimising weak lensing survey 
designs. For instance, switching from the COSMOS non- 
LRG sample (which is typical of a 'normal' weak lensing 
sample) to one similar to the LRG sample would yield the 
same shape noise level if the number density of the latter 
were about a factor of 2 smaller (ignoring the contribution by 
measurement noise which can be independently controlled, 
e.g. by apparent magnitude cuts). 

The actual LRG sample used in this work contains only 
a small percentage of the total number of galaxies suitable 
for weak lensing measurement, but the Iow-CTc found also in 
part of the non-LRG sample indicates that there is scope for 
a more optimal cut that maximises the difference in ae while 
balancing the sample sizes. However, the choice of an opti- 
mal galaxy sample for the cosmological exploitation of weak 
lensing is more complex as, in addition to the intrinsic ellip- 
ticity dispersion, it needs to take into account the quality of 
photometric redshifts attainable, the anticipated constraints 
on intrinsic alignment contamination, the presence of colour 
gradients across galaxy images, etc. 

In spite of the necessarily still simplistic assumptions 
about the shapes of galaxies and how these are linked to 
the dark matter properties, the best-faring simulation-based 
models were capable of yielding fair agreement with the ob- 



served (Te . The difference in the amplitude of ae between the 
LRG and non-LRG sample is reproduced quantitatively, in- 
cluding the independence on redshift. The models do not 
show the drop in ae for luminous late-type galaxies, which 
however might be accounted for by allowing a variation of 
the ratio of disk thickness to size r edge-on with galaxy pr op- 
erties (e.g. to match the results of lBailin fc Harrisll2008r ). 

While producing distributions of |e| that are reasonably 
close to the observations for the LRG sample, all models 
yield discrepant results for the corresponding distribution 
of the non-LRG sample, particularly at low ellipticities. As 
discussed above, introducing a distribution of non-circular 
light profiles for disks might resolve this issue. Observa- 
tions, which should be feasible at least in the local Uni- 
verse, could inform future implementations of models about 
both the deviations from circular disks and the distribution 
of rcdge-on . Incorporating the shape properties of irregular 
galaxies, which constitute a substantial fraction of the full 
sample at redshifts around unity, is however an open issue. 

The comparison between models and COSMOS data 
suggests that a small value of r-cdgo-on, i.e. a subdominant 
bulge, is preferred. However, this may at least in part be 
an artifact of the failure to produce the correct fraction of 
low-ellipticity galaxies in the models. Furthermore it is evi- 
dent that early-type galaxy models based on reduced inertia 
tensor measurements of halo shapes do not match the obser- 
vations at alo. This is a surprising finding, given that the 
reduced inertia tensor assigns more weight to the interior of 
the halo, i.e. to regions closer to where the galaxy resides. 

The approach taken in this paper can readily be applied 
to higher resolution simulations, where shapes and angular 
momenta of individual haloes are obtained also for satel- 
lite galaxies. Further improvements are expected from more 
observations of ensemble properties of galaxy shapes (such 
as for rcdgc-on, see above), as well as from more and larger 
hydrodynamic simulations which give insight into the con- 
nection between the shapes of the dark matter halo and the 
luminous galaxy. 

But even at the present stage, the analysis of one-point 
statistics of galaxy shapes in the form of the polarisation dis- 
persion has proven to be insightful and capable of discrimi- 
nating between models with high significance. Compared to 
measurements of correlations of intrinsic galaxy shapes, dis- 
persions and distributions of e can be obtained deeply into 
high-density (non-linear) regions, are not affected by galaxy 
bias, and are less susceptible to sample variance. Hence, to 
constrain models of intrinsic galaxy shapes and alignments, 
the approach presented here is a valuable complement to 
second-order statistics of galaxy shapes. 
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Note that our procedure is not equivalent to simply pro- 
jecting the inertia tensor M along the line of sight as this 
would only correspond to the first term contributing to W~^ 
in Equation (|A4[) . In other words, the shape derived from 
the inertia tensor of the two-dimensional mass distribution 
on the sky is generally not the same as the projected shape 
of the ellipsoid derived from the three-dimensional inertia 
tensor. This ambiguity arises because of our simplistic as- 
sumption of a top-hat radial light distribution, in one case 
implicitly imposed on the two-dimensional, and in the other 
on the three-dimensional galaxy model. We choose to imple- 
ment the latter approach as it is more versatile with respect 
to future implementations of the galaxy shape model, e.g. 
the introduction of a misalignment between the major axes 
of dark matter halo and galaxy. 

Using the definitions given in Section [3. II the eigende- 
composition of M can be written as 

M = VD V", where T:",,^ = S^puj^-Y = {3^32,83}. (A5) 

Inserting these expressions into Equation (|A4|) . and identi- 
fying a = ri^ M~^ n as well as fe = M~^ n, it is straight- 
forward to derive Equations Q and Q. 



APPENDIX A: PARALLEL PROJECTION OF 
AN ELLIPSOIDAL HALO 

We define the shape of a dark matter halo and the corre- 
sponding early-type galaxy in terms of the eigenvectors and 
eigenvalues of the halo inertia tensor M, given by Equation 
lU. Therefore the surface of the ellipsoid that serves as the 
model of the gala:xy is constituted by all points with coor- 
dinates X which fulfil 



c^ M"^ X = 1 



(Al) 



The line of sight is chosen to coincide with the third coor- 
dinate axis along which we parallel project the ellipsoid. To 
this end we define sets of straight lines 

x{t) = y + tn with n = {0,0,1}"; y = {yl,0}", (A2) 

i.e. these lines are parallel to the line of sight and, for t = Q, 
intercept the plane onto which the ellipsoid is projected at 
position y^. 

Inserting Equation (|A2|) into Equation (|A1[) . one ob- 
tains a quadratic equation in t of the form 



Ar + Bt + C = with 



(A3) 



A = n MT'^ n ; B = 2nM.'^y- C ^ y^ M"^ 3/ - 1 • 

Those lines that touch the surface of the ellipsoid in one 
point are selected by requiring that the discriminant van- 
ishes, 4AC — B^. This condition is equivalent to 



y^M-^y 



= y 



y 



M 



w 



(n-M-'y)- 



(A4) 



yi 



In the last line we have defined the symmetric two- 
dimensional tensor W. This expression is the defining equa- 
tion of an ellipse that lies in the plane of the sky. 



APPENDIX B: EXTRACTING THE INTRINSIC 
POLARISATION DISPERSION FROM WEAK 
LENSING SHEAR CATALOGUES 

Bl Fisher measurement noise correction 

To estimate the contribution by measurement noise to the 
observed polarisation dispersion, we follow the method pro- 
posed bv lLeauthaud et al.l ([20071). Galaxy images are mod- 
elled as a bivariate Gaussian, 



G{x) 



S 



, exp < —-(x — mV Q ^ (x — m) > , (Bl) 

27r^/ditQ ^l 2^ J ^ y yj-A ; 

where S denotes the fiux, m the position of the image centre, 
and Q the symmetric second-order brightness tensor. These 
six parameters will in the following be collected into a six- 
dimensional parameter vector p. 

The model G{x) can be fit to the actual light distri- 
bution of a galaxy image. To avoid working at the image 
level, one resorts to the Fisher matrix which is given by 
the expectation value of the Hessian of the log-likelihood 
corresponding to this fit. This Fisher matrix is given by 
(|Leauthaud et al.ll2007l ') 






2 / 



dG(x,) dGjXi) 



(B2) 



where (tjv is the noise per image pixel. The sum runs over all 
pixels covered by the image, which in our implementation 
comprises the pixels within twice the half-light radius. Note 
that the derivatives in Equation (|B2|) are readily computed 
analytically. 

Assuming that the image centre is at m = 0, the only 
ingredients needed to compute the Fisher matrix are the 
flux S, the noise level (Tjv which can be inferred from the 
SExtractor FLUXERR_AUTO parameter, and the three 
components of the brightness tensor Q. The shear catalogue 
that we are working with provides us only with the polari- 
sation components 61,2, i.e. the image size information con- 
tained in Q is not directly available. 
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We retrieve this information in an approximate fashion 
by calculating the quantity 

_ Qii + Q22 _ J d^a; Kksb{x) I{x) \x\^ 



- r , (B3) 

2 2 J d^x Kksb{x) I{x) 

where -R'ksb is the KSB weighting, which ISchrabback et alj 
l|201Cl ) chose as a circular Gaussian with the half-light radius 
as width. Note that the introduction of a weight function in 
the brightness moments is essential, as otherwise Q would be 
dominated by noise in the outskirts of the image. We model 
the galaxy's light distribution I{x) as a circular Sersic profile 
with a typical Sersic index between 1 and 4 and a scale radius 
given in terms of the half-light radius. Then the elements of 
the brightness tensor are given by 



= s(l 



ei 



>i2 = se2 ; Q22 = s(l 



ei 



(B4) 



The flux S required by the model in Equation (IB1|I is 
not necessarily equal to the measured flux of the galaxy 
image, as it is the best-fit amplitude of the Gaussian 
model. These two quantities can differ substantially since 
the Gaussian model do es not provide very accurate fits. 
iLeauthaud et al.l (|2007l ) introduced an overall calibration 
factor to account for this and determined it via image simu- 
lations. Instead, we choose an analytical route and calculate 
a calibration factor for each galaxy individually as the best- 
fit amplitude of a radial Gaussian profile to a radial Sersic 
profile, multiplied by the KSB kernel (both depending on 
the half-light radius of the galaxy). Note that this ansatz 
again makes the simplifying assumption of circular galaxy 
images. 

After performing these steps, the Fisher matrix is com- 
puted for each galaxy. The submatrix corresponding to the 
three Q elements of the inverse Fisher matrix yields an es- 
timate of the covariance Cov(Q) of Qii, Q12, and Q22- E)ue 
to taking the inverse of the full six-dimensional Fisher ma- 
trix, the covariance is marginalised over uncertainties in the 
centroid position and the flux. The measurement error on 
lel^ can then be derived via 



Alel- 



\ 



d\ef 



3 

^ 'dQl 



[Cov(g)], 



dQi 



(B5) 



where we introduced combined indices i,j G {11; 12; 22}. 
The derivatives are obtained in analytic form by making use 
of the d efinition of the polarisation in t erms of the brightness 
tensor (|Bartelmann fc Schneideill200H ). The error according 
to Equation (|B5|I is computed for every galaxy, averaged 
over the samples and bins under consideration, and sub- 
tracted from the dispersion of e obtained via Equation (|10|) . 
Our computation of the brightness tensor elements re- 
quires |e| ^ 1, but in practice this limit can be exceeded for 
estimators of e due to noise. This happens for 84 galaxies 
in our sample (all but one with _F814M^ > 24), which we 
discard completely. 



B2 Circularisation correction of the polarisation 

To determine the effect of the circular Gaussian kernel in- 
cluded in the brightness moment s in the KSB implemen- 
tation of ISchrabback et al.l l|2010r ) on the measurement of 




Figure Bl. Relation between the absolute value of the polarisa- 
tion as measured by KSB, Cobsi ^ind the true polarisation etruo- 
The blue (red) line results when assuming a galaxy light profile 
with Sersic index nscrsic = 1(4)- The black dotted line indicates 
a one-to-one relation. 



galaxy ellipticity, we resort again to analytic light distribu- 
tions using Sersic profiles. Varying the image polarisation 
between zero and unity (assuming all isophotes have the 
same polarisation), we compute the circular half-light ra- 
dius for each profile. The result is then fed into the Gaussian 
kernel A'ksb, which, together with the light distribution, is 
used to calculate the sec ond brightness moments. By rn eans 
of the defining equation l|Bartelmann fc Schneiderll2001f) the 
'observed' polarisation is derived from the brightness tensor. 

In Fig. lBll the input (and hence true) absolute value of 
the polarisation in our computation, termed et,uc, is plotted 
against the resulting eobs for a Sersic index of nscrsic = 1, 
appropriate for late-type galaxies, and nscrsic = 4, corre- 
sponding to a de Vaucouleurs profile typical of early-type 
galaxies. The observed polarisation is considerably smaller 
than the true galaxy polarisation, caused by the circular 
weighting entering the brightness moments. The dependence 
on the Sersic index is weak, so that we correct all galaxies 
according to a curve with nsersic = 3.25, which is close to 
the average of the two lines shown in the plot. 

In our alternative approach to determine Oe the noise- 
corrected polarisations are subjected to the correction for 
the Gaussian kernel circularisation, based on the relation 
shown in Fig. lBll This shifts the violet dotted line of Fig.|4] 
up to the black dotted line by about 0.15 in CTe, indepen- 
dent of apparent magnitude by construction. Note that this 
correction is included in the 'shear tensor' of the KSB for- 
malism and thus automatically accounted for in the shear 
estimates which form the basis of our default approach. 

Such a large correction may be regarded as a strong ar- 
gunrent against employing weak lensing shear estimates for 
intrinsic galaxy shape measurements. However, it should be 
kept in mind that re-using these catalogues is not only con- 
venient as all necessary steps to eliminate PSF effects and 
other systematics have already been performed, but weak 
lensing shape measurement methods might well be the only 
way to obtain reliable intrinsic shapes in the low signal-to- 
noise and small apparent s ize regime. Wheth er approaches 
other than KSB (see e.g. iBridle et al.l l2010f ) are perhaps 
more suitable for this purpose remains the scope of future 
work. 
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